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Abstract. In [3], the authors proved an existence result for BSDEs with quadratic 
generators with respect to the variable z and with unbounded terminal conditions. 
However, no uniqueness result was stated in that work. The main goal of this paper 
is to fill this gap. In order to obtain a comparison theorem for this kind of BSDEs, 
we assume that the generator is convex with respect to the variable z. Under this 
assumption of convexity, we are also able to prove a stability result in the spirit of 
the a priori estimates stated in [6]. With these tools in hands, we can derive the 
nonlinear Feynman-Kac formula in this context. 

1. Introduction. 

Since the first existence and uniqueness result for backward stochastic differential equations 
(BSDEs in short for the remaining of the paper) of Pardoux and Peng [13], lots of works have 
been done in this area and the original Lipschitz assumption on the generator, the function / 
in the BSDE 

Y t =£ + J f(s,Y s> Z s )ds- J Z s -dB s , 0<t<T, (1) 

has been weakened in many situations. Let us recall that, in the previous equation, we are 
looking for a pair of processes (Y, Z) which is required to be adapted with respect to the 
filtration generated by the Brownian motion B. Even though there exist some improvements 
of Pardoux-Peng's result for multidimensional BSDEs (see e.g. [12] or [2]), the case of real 
valued BSDEs meaning that Y is a real process is easier to handle. The main reason for the 
gap between the real case and the multidimensional one comes from the comparison theorem 
which turns out to be a powerful tool to construct solutions for real valued BSDEs. Roughly 
speaking, a comparison theorem for BSDEs allows one to compare solutions to BSDEs as soon 
as one can compare the terminal conditions (£ in the previous equation) and the generators. 
In particular, a classical strategy in order to obtain a solution to the BSDE (1) consists in 
constructing an increasing sequence of solutions to BSDEs and then passing to the limit. This 
method is used for instance in [10] for continuous generators with a linear growth. But one of 
the most important application of this approach is the work of M. Kobylanski [9] concerning 
quadratic BSDEs. We should point out that quadratic BSDE means a BSDE whose generator 
has a quadratic growth with respect to the variable z. For these quadratic BSDEs, all the 
classical results, existence and uniqueness, comparison theorem and stability of solutions, has 



1 



last modif. 2007-03-14 



been stated in [9] but with the restriction that the terminal conditions have to be bounded 
random variables. 

It was known for a long time that the boundedness of the terminal condition is not really 
needed for studying a quadratic BSDE. Actually, on simple examples, on can see that the 
existence of an exponential moment of sufficiently large order is enough to get a solution. In 
[3], existence of solutions to quadratic BSDEs with unbounded terminal conditions was studied 
and it was proved that the existence of an exponential moment of £ is a sufficient condition in 
order to construct a solution. However the question of uniqueness of solutions in this setting was 
left open in that work. Quadratic BSDEs with bounded terminal conditions has an interesting 
feature : Z • B, the stochastic integral of the process Z with respect to the Brownian motion 
B, is a BMO-martingale. This observation can be used to obtain uniqueness see [8] and also 
other properties [1]. When the terminal condition is not bounded, Z*B is, in general, no more 
a BMO-martingale: this is a difficulty for uniqueness. 

The main contribution of this paper is to obtain a uniqueness result for quadratic BSDEs 
with unbounded terminal conditions when the generator of (1) is assumed to be in addition 
convex or concave with respect to the variable z. Let us mention that convex generators appear 
naturally for BSDEs associated to stochastic control problems, see e.g. [7]. With this further 
assumption, it is possible to obtain a comparison theorem for solutions to unbounded quadratic 
BSDEs which of course gives uniqueness see Section 3. The key idea for proving this result is 
the following: instead of trying to estimate the difference of two solutions, say Y and Y', we 
estimate, for each £ (0, 1), Y — BY' \ this allows to take advantage of the convexity of the 
generator. Moreover, it turns out that the convexity of the generator with respect to the variable 
z is also a convenient assumption to obtain a stability result. The proof of this last result relics 
mainly on the same computation mentioned before even though technical difficulties arise and 
impose us to go into details carefully. Finally all these results yield the nonlinear Feynman-Kac 
formula in this framework. 

The paper is organized as follows. In Section 2, we prove an existence result in the spirit of 
[3]: here, we work with generators which are quadratic with respect to z and monotone with 
respect to y. This monotonicity assumption allows as usual to get rid of the growth assumption 
on the generator with respect to y. Section 3 is devoted to the comparison theorem from which 
we get as a byproduct an existence and uniqueness result for quadratic BSDEs with unbounded 
terminal conditions. In Section 4, we obtain a stability result in our framework and finally in 
the last section we derive the nonlinear Feynman-Kac formula. 

Let us close this introduction by giving the notations which we will use in all the paper. For 
the remaining of the paper, let us fix a nonnegative real number T > 0. First of all, B = {B t }t>o 
is a standard Brownian motion with values in R d defined on some complete probability space 
(CI, J 7 , P). {Ft}t>o is the natural filtration of the Brownian motion B augmented by the P-null 
sets of J- '. All the measurability notions will refer to this filtration. In particular, the sigma-field 
of predictable subsets of [0, T] x fl is denoted V . 

As mentioned in the introduction, we will deal only with real valued BSDEs which is an 
equation of type (1). The function / is called the generator and £ the terminal condition. Let us 
recall that a generator is a random function / : [0, T] x O x R x R d — > R which is measurable 
with respect to V <E> B(R) <%) B(R d ) and a terminal condition is simply a real .^-measurable 
random variable. 

By a solution to the BSDE (1) we mean a pair (Y, Z) = {(Y t , Z t )}t£[o,T] of predictable 
processes with values in R x R d such that P-a.s., t i — > Y t is continuous, t i — > Z t belongs to 
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L 2 (0,T), 1 1 — ► f(t,Y t ,Z t ) belongs to 1/(0, T) and P-a.s. 

Y t=Z + J t f{s,Y s ,Z s )ds- J Z s -dB s , 0<t<T. 

We will sometimes use the notation BSDE(£, /) to say that we consider the BSDE whose 
generator is / and whose terminal condition is £; (V^(£), Z?((;)) means a solution to the 
BSDE(£, /). A solution (Y f (£), Z^ {£)) is said to be minimal if P-a.s., for each t e [0,T], 

Y/(£) < Yf(C) whenever P-a.s. £ < C and f(t,y,z) <g(t,y,z) for all (t,y,z). (^(0,^(0) 
is said to be minimal in some space B if it belongs to this space and the previous property holds 
true as soon as (F ff (C), Z 9 (C)) G £. 

For any realp > 1, <S P denotes the set of real- valued, adapted and cadlag processes {^t}te[o,T] 
such that 

||y|| 5P :=E[sup < t < T |y t pf /p <+oo. 

(S p , || ■ H^p) is a Banach space. 

M p (R d ) or simply M p denotes the set of (equivalent classes of) predictable processes 
{Zt}te[o.T] with values in R d such that 

\\Z\\ M p :=e 

M p (R rf ) is a Banach space endowed with this norm. We set S = U p> iS p , M = U p >iM p and 
5°° stands for the set of predictable bounded processes. 

Finally, we will say that a real process Y, adapted and cadlag, belongs to £ if the random 
variable Y* = sup 0<t<T |Yt| has exponential moments of all orders and we recall that Y belongs 
to the class (D) as soon as the family {Y T : t < T stopping time} is uniformly integrable. 

2. An existence result. 

In this section, we prove a mere generalization of the existence result for quadratic BSDEs we 
obtained in [3] . We consider here the case where the generator has some monotonicity property 
with respect to the variable y. As usual, this kind of assumptions allows to get rid of the (linear) 
growth of the generator / in y. See e.g. [12], [2] or [4]. Of course, we still consider the case of a 
generator with a quadratic growth in the variable z and of an unbounded terminal condition. 
In this section we assume the following on the generator. 

Assumption (A.l). There exist two constants (3 > and 7 > together with a progres- 
sively measurable nonncgative stochastic process {a(i)}o<t<T and a deterministic continuous 
nondecreasing function ip : R + — > R + with </?(0) = such that, P-a.s., 

(i) for all t £ [0, T], (y, z) \ — > f(t, y, z) is continuous; 

(ii) monotonicity in y : for each (t, z) G [0, T] x R d , 

VyGR, y(f(t,y,z)-f(t,0,z))<p\y\ 2 ; 

(iii) growth condition: 

V(f, y, z) G [0, T] x R x R d , \f(t, y,z)\< a(t) + p(|y|) + l\z\ 2 . 
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2 )dt 



Let us denote by h the (random) function h(t, x) = a(t) + (3x. First of all, we want to 
derive a sharp estimate for solutions to the BSDE (1). The idea is the following: we apply 
Ito-Tanaka's formula to compute U(t, \Y t \) = e r ^^ Y *^' . Let us denote L the local time of Y at 
0. We have 

\Y t \ = \£\+ [ sgn(Y s )f(s,Y s ,Z s )ds~ [ sgn(Y s )Z s ■ dB s - [ dL s , 
Jt Jt Jt 

and, from the growth of / (A.l)(ii)-(iii), 

sgn(F s )/( S , Y„ Z s ) = sgn(r s ) (f(s, Y s , Z s ) - f(s, 0, Z s j) + sgn(Y s )f(s, 0, Z s ) 

< 0\Y.\ + a(s) + 1\Z S \ 2 = h(s, \Y„\) + \\Z S \\ (2) 

From Leo's formula, we derive the inequality 

f^f]} = {-Mt, \Y t \) sgn(Y t )f(t, Y t , Z t ) + Ml \Yt\) + |tfx(<, \Yt\f\Z t 

+ ^xx(t, \Y t \)\Z t \ 2 dt + 4> x (t, \Y t \)dL t + ^ x {t, \Y t \)sgn(Y t )Z t ■ dB t . 
If ipx(t,x) > 1 for x > 0, we have, taking into account the inequality (2), 
Mt, \Y t \)sgn(Y t )f(t,Y t , Z t ) - Mt, \Yt\) ~ \^ \Yt\f\Z t \ 2 < r/, x (fi, \Y t \)h(t, \Y t \) - <fc(t, \Y t \). 

The idea is now clear; we want to find a solution ip to the PDE: for s E [0, T], 

i>t{t,x)-h(t,x)^ x (t,x) =0, ip{s,x) = *(x), tE[s,T], 

such that ijjx {t, x) > 1 and ip X x(t, x) > for x > 0. But the previous PDE can be solved by the 
characteristics method. Let {v(u;t, x)}o< u <t be the solution to the integral solution 

v(u\ t, x) = x + / h(r,v(r;t,x)) dr, 0<u<t. 
J u 

Then i/}(t,x) = \& (i>(s; t, x)). In our case, h{r,x) = a{r) + [3x, the solution is given by the 
formula 



v(u; t, x) = xe p{t - u) + / a{r)e 0{r - u) dr. 

J U 

We choose ^f(x) = x and we get, if a is a stopping time such that s < a < T, 

e^ Y '\ = U(s, \Y S \) < U(a, \Y a \) - f jU(t, \Y t \)^ x (t, \Y t \) sgn(Y t )Z t ■ dB t . (3) 



This computation leads the following a priori estimate. 
Proposition 1. Let (Y,Z) be a solution to the BSDE (1) such that the process 



exp [ ~ie l3T \Y t \ + 1 J a(r)e fjr dr j 
belongs to the class (D). Then, P-a.s., 

VtE[0,T], \Y t \ < ilogE ^exp ^ 7e / s C T -*)|f| + 7 ^ a(r)e' 3(r -* ) dr^ T t ) . (4) 
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Proof. Let s E [0,T]. Set, as before ip(t,x) = xe^'^ + f g a(r)e /3 ( r " s ) dr, U{t,x) = e^*'*). 
Let us consider, for each integer n > 1, the following stopping time: 

o-n = inf jt G [s,T] : 7 2 ^ £/ 2 (r, \Y r \)i/) x (r, \Y r \f\Z r \ 2 dr > n| A T. 

It follows from the inequality (3) and the definition of <t„ that 

e 7ini < E (>v(^,i^j) i ^ < E |^ cxp / 7 ^(t- s) | Fct j + 1 J\ {r y^s )d ^ | ^ _ 

Thus, the inequality follows by sending n to infinity. □ 

Remark. It is clear, from the previous computations, that an a priori estimate can be stated 
for functions h which are not linear as we did. For instance, it is possible to obtain such an 
estimate when 

y(f(t,y,z)-f(t,0,z))<\y\p(\y\) 
where p : R + — ► R + is convex, C 1 with p(0) = 0; \ f(t,y,z)\ < a + (p(\y\) + §|z| 2 (a > 0) and 

/ — — - = +oo, with h(u) = p{u) + a. 
Jo h(u) 

In this case, one has to choose i()(t, x) = _1 (t — s + 0(x)) where 

8(x) 

See e.g. [11] or [3]. 



h(u) 



As a byproduct of this a priori estimate, we see that, when |£| and |a|i := J Q a(s)ds are 
bounded random variables, any solution (Y, Z) to (1) such that e ie<i y * belongs to the class 
(D), is actually bounded. More precisely, 

im^mei + Miiioo. 

This observation allows to generalize a little bit, with exactly the same proofs, some of Koby- 
lanski's results [9]. 

Lemma 2 (M. Kobylanski). Let the assumption (A.l) hold. If ^ and \ot\i are bounded random 
variables, then the BSDE (1) has a minimal bounded solution (and a maximal one also). 

Moreover, let {f n )n>i be a sequence of generators satisfying (A.l) with (a, /3, / y, <p) indepen- 
dent of n such that (f n (t, y n , z n)) n >i converges to f(t,y,z) as soon as (y n ,z n ) — ► (y, z) and 
1st (£n)n>i be a sequence of terminal conditions converging almost surely to £. Let us assume 
that, for each n > 1, the BSDE associated to and f n has a solution (Y",Z n ) G <S°° x M 2 
such that (y n )n>l * s nondecreasing and, for some constant C, sup ra>1 ||F ra ||oo < C. 

Then, there exists (Y, Z) G 5°° xM 2 , solution to the BSDE (1), such that (Y n ) n >i converges 
to Y uniformly on [0,T] in probability and (Z n ) n >i converges to Z in M 2 . 

We are now in position to generalize to this framework the result we obtained in [3] . 
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Proposition 3. Let (A.l) hold. 

If I CI + Mi has an exponential moment of order 7e^ T , then the BSDE (1) has a solution 
(Y, Z) such that 

Vt€[0,T\, \Y t \ < ilogE ^exp ^^"'^l +1 a{r)e^ r -^dr\ | ^ J . (5) 

Proof. Let us assume that £ and / are nonnegative. For each integer n > 1, we consider 
£„ = £An and f n (t,y,z) = lt<o n f{t,y,z) where 



cr„ = inf jf e [0, T] : J a{s) ds > n| A T. 



According to the first part of the Lemma 2, let (Y n , Z n ) <E 5°° x M 2 be the minimal bounded 
solution to the BSDE whose terminal condition is £„ and whose generator is f n . It follows, 
from Proposition 1, that for each n > 1, 

\Y t n \ < 1 log! (cxp (7^ T + Mi)) | T t ) := X t . 

Since £ n < /n £ /n+i an( l since we are dealing with minimal solutions, the sequence (Y n ) 

is nondecreasing. We have in mind to use the stability property quoted in the second part of 
Lemma 2. To use this result, we need to know that the sequence (Y n ) remains bounded which, 
of course, is not true in general. To overcome this difficulty the idea is to use the previous 
estimate and to work on random time interval instead of working on the whole interval [0, T]. 
To do this, let, for k > 1, Tfe be the following stopping time 

r fc = mf{t G [0, T] : X t > k} A T. 

By construction the sequence (in n) Y k (t) = Y™ hT remains bounded in k. Setting moreover 
= Z?l t < Tk , we have 

Y k n (t) = F™ + ^ l t < Tk ^J (s, Y k n (s), Z%(s)) ds - ^ Zl{s) ■ dB s . 

We apply the second part of Lemma 2 to obtain, for each k, a solution (Y k , Z k ) to the BSDE 

Y k (t) = a + j " f (s, Y k {s),Z k {s)) ds- Z k (s) ■ dB s , with & = sup^ Y T n k . 

But by construction, r k < r k +\ so we have the localization property 

Y k+1 (t A r k ) = Y k {t) 7 Z k+1 {t)l t < Tk = Z k {t). 
If we define the processes Y and Z on [0, T] by the formula 

Y t =Y 1 (t)l t < Tl +J2 Y k(t) 1 ]T k - 1 ,r k ](t), Z t = Z 1 {t)l t < ri +Y l Z lc (t)l ]Tk _ uTk] {t), 

k>2 k>2 

the previous BSDE is rewritten as follows 



Y t = b, + J " f(s, Y„ Z s ) ds- J" Z s - dB s . 
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The last point is the fact that, P-a.s., Tfc = T for k large enough which allows to send k to 
+00 in the previous equation to prove that (Y, Z) is a solution to (1). Of course, the inequality 
of the theorem is verified by the process Y since it is verified by each process Y n in view of 
Proposition 1. 

In the general case, we use a double approximation; £™' p = £ + A n — £~ A p, 
f n '"(t,y,z) = l t <aJ + (t,y,z)-l t < !7 J-(t,y,z). 

□ 

Corollary 4. Let the assumption (A.l) /10/d. 

// |£| + |a|i /ias an exponential moment of order A > 7e^ T t/iera t/ie BSDE (1) /ias a solution 
(Y,Z) such that Y satisfies the inequality (5) and Z belongs to M 2 . 

In the case where |£| + |a|i has exponential moments of all order, this solution is such that 
Y belongs to £ and Z belongs to MP for all p > 1. More precisely, for each p > 1, 

exp Ctp sup \Y t \] + ( f \Z s \ 2 ds\ 

\ 0<t<T / \Jo I 
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<<7E[exp(p 7 (|£| + |a|i))] (6) 



where C depends on p, 7, (3 and T . 



Proof. Let (Y, Z) the solution to (1) obtained in Theorem 3. Then Y satisfies the estimate (5). 
It follows immediately from this estimate and Doob's maximal inequality that e lY belongs to 
S q for some q > 1 when £ := |£| + |a|i has an exponential moment of order A > 7e^ T and that 
the estimate (6) holds true for Y when £ has exponential moments of all order. To obtain the 
results for the process Z, we argue as in [3]. 

For n > 1, let r„ be the following stopping time 

inf<fi>0: [ e 2 ^\Z s \ 2 ds > n\ AT, 



and let us consider the function from R + into itself defined by 

u(x) = [e lx — 1 - 7x) . 
7 

x 1 — > u(\x\) is C 2 and we have from Ito's formula, with the notation sgn(x) = —l x <o + lx>0: 

u(\Y \) = u (|y tAT j) + " ^'(|y s |)sgn(y s )/( s ,y s ,z s )-i u "(|y;|)|z s | 2 ^ s 

*At„ 

U '(|y s |)sgn(K,.)Z s -rfB.,. 
It follows from (A.l) that 

sgn(F s )/( S , Y s , Z s ) < sgn(r s ) (f(s, Y s , Z s ) - f(s, 0, Z s )) + \f(s, 0, Z s )\ < a(s) + (3\Y S \ + 1\Z S \ 2 , 
and since u'(x) > for x > that 

u{\Y \) <u(\Y tATn \)+ I u'{\Y a \){a(s)+f3\Y s \)ds- \ u'(\Y s \) sgn(Y s )Z s ■ dB s 
1 

~ 2 



(At„ 

K(|F s |)-7 U '(|F s |))|Z s | 2 d S . 



r 



Moreover, we have (u" — ju')(x) = 1 and u(x) > for x > 0, so taking the expectation of the 
previous inequality 







f " \Z s \ 2 ds 


< E 


Jo 





-~ sup e 

7 te[o,T] 



7ini , ± 



7 Jo 



(a(s) + p\Y a \) da 



Fatou's lemma together with the fact that e 7 l y * G S q for some q > 1 gives the result when £ 
has an exponential moment of order A > r ye" T . 

For the second part of the result, let us observe that 



1 



Z s \ 2 ds<— supe 7|yt| + -/ e^ y ^ (a{s) + f3\Y s \) ds 



e i\Yt\ 

7 Jo 

r*AT„ 

u '(|r s |)sgn(y s )z s -dB £ 



sup 

0<t<T 



II 



It follows from the BDG inequalities, since \y\ < e 7 ' J, '/7 an d U '(M) ^ e 1 ^/^, that 



E 



p/2 



< CE 



< C'l 



& PlY* +eP7 y/2| a |P/2 + C / " e 2 7 |Y s || Zs| 2 ds 



P7^* +e f 7Y */2| Q |P/2 



LZ S r/.s 



P/4' 



p/2 



from which the result follows using Fatou's lemma. 



□ 



3. Comparison Theorem. 

Let us consider now the main topic of this paper : uniqueness for quadratic BSDEs with 
unbounded terminal condition. We state in this section a comparison theorem for solutions to 
quadratic BSDEs and let us recall first that, in the bounded case, such a result exists (see [9]) 
even though it requires more assumptions than those in the existence result. In our unbounded 
framework the situation is exactly the same : the assumptions to prove uniqueness arc stronger 
than those needed for existence. In particular, a very convenient way to derive the comparison 
theorem is to assume that the generator is convex (or concave) with respect to the variable z. 
From the point of view of PDEs, the convexity of the generator seems a natural assumption see 
e.g. [5]. Let us consider the following assumption on the generator /. 

Assumption (A. 2). There exist two constants 7 > and (3 > together with a nonnegative 
progressively measurable stochastic process {ct(t)}o<t<T such that, P-a.s., 

(i) for all t £ [0, T], for all y € R, z 1 — ► f(t, y, z) is convex; 

(ii) for all (t, z) G [0, T] x R d , 

V(w')eR 2 , \f(t,y,z)-f(t,y',z)\ <p\y-y'\; 

(iii) / has the following growth: 

y, z) G [0, T] x R x R d , \f(t,y,z)\ < a(t) + /3 \y\ + ^ \z\ 2 ; 

(iv) |a|i has exponential moment of all order. 



<s 



Theorem 5 (Comparison theorem). Let (Y, Z) be a solution to (1) and (Y',Z') be a solution 
to the BSDE associated to the terminal condition £' and to the generator f such that both Y 
and Y' belongs to £, namely, 



e XY* + e XY" 



< + OO. 



VA > 0, E 

We assume that, W-a.s., 

t<?, y(t, y, z) e [0, T] x R x R d , f(t,y,z)<f%y,z). 

If f verifies (A.2), then F-a.s., for each t e [0,T], Y t < Y{. 
If moreover, Yq = Yq, then 

(f'-f)(t,Yj,Z' t )dt = 0) >0. 



Proof. The idea is the following: instead of trying to estimate the difference between the pro- 
cesses Y and Y' , we estimate Y — 9Y' , for each 9 £ (0, 1), in order to take advantage of the 
convexity of the generator. Similar idea is also used for quadratic PDEs in [5] where the authors 
prove that u — 9u' is nonpositive for two solutions u and u' . 

Let 9 € (0, 1) and let us set Ut = Y t — 6Y( and Vt = Zt — 9Z[. Let us consider a real 
stochastic process {a(t)}o<t<T, progressively measurable, with integrable paths to be chosen 
later. We set, for all t E [0, T], A t = f Q a(s) ds. Wc have, from Ito's formula, 

e At U t = e AT U T + J e A °F s ds - J e As V s ■ dB Sl 0<t<T, 
where, of course, 

F t = (f(t, Y u Z t ) - 9f (t, Yl, Z' t )) - a(t)U t . (7) 
Let us set as usual Sf(t) = (/ — /') (t, Y(, Z' t ) so that 

f(t, Y t , Z t ) - Of (t, Yl, Z't) = f{t, Y t , Z t ) - Of it, Yl, Z' t ) + 95 fit) 

= (f{t, Y t , Z t ) - / it, Yl, Z t )) + (f (t, Yl, Z t ) - 9f (t, Yl, Z' t )) + 95 fit). 

Since / is convex with respect to z, the second term of the right hand side of the previous 
inequality can be easily bounded from above. Indeed, 



fit,Yl,Z t ) = f[tX,9Z' t + i\ 



' 1-9 

< 9f it, Yl, Z' t ) + (1 - 9)f (t, Yl, 

and from the growth of the generator /, (A.2)(iii), 

/ it, Yl, Z t ) < Of it, Yl, Z't) + (l-9) (ait) + (3 + \Z t - OZ'f . (8) 

Roughly speaking, the first term does not change anything since it depends more or less 
only on y and / is Lipschitz with respect to this variable. We get rid of it with a classical 
linearization. Let us write 

/ it, Y t , Z t ) - f it, Y t ', Z t ) = f it, Y t , Z t ) - f (t, OYl, Z t ) + f (t, OYl, Z t ) - f (t, F/, Z t ) 

= a(t)U t + f (t, OYl, Z t ) - f (t, Yl, Z t ) 
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where a(t) = [/ (i, Y t , Z t ) — f (t, 6Y(, Z t )] jXJt when Ut ^ and a(t) = (3 in the other case. Since 
/ is /3-Lipschitz, a is bounded by (3 and 

/ (t, Y t , Z t ) - / (t, Y/, Z t ) < a(t)U t + (1 - 6)p |Y/| . (9) 

If we choose for a the process we have just introduced, we get, from (8) and (9), coming 
back to the definition of F (7), 

Ft < (1 - 9) (a(t) + 2/3 |Y/|) + ^L— \V t \ 2 + 95f(t). (10) 

Now, we get rid of the quadratic term with an exponential change of variables. Let c > 
and let us set Pt = e ce tUt , Qt = cP t Vte At . From Ito's formula we deduce that 

Pt=P T + c£ P s e A ^F s -^\V s \ 2 ^jds-c^ e A °P s V s -dB s 
:=P T + [ G s ds- [ Q s -dB s . 



Jt Jt 

(10) yields, since c is nonnegative, 



G t < cP t e A * {(1 - 9) (a(t) + 2(3 |Y/|) + 95f(t)} + cP t e A * (^TJZ^) ~ M 

But A t > -j3T so that, if we choose c = c(9) := 7e /3T /(l - 9) we obtain the following inequality, 
G t < P t e At (6c(0)6f(t) + ~fe pT (a(t) + 2(3 |Y/|)) . (11) 
Finally, let us introduce the processes 

D t = cxp Qf e A ° (9c(9)Sf(s) + 7 e^ T (a(s) + 2/3 |Y S '|)) ds\ , P t = D t P t , Q t = D t Q t . 
Once again Ito's formula gives us, for any stopping time r such that < t < r < T, 

Pt<Pr- I Qs-dB s . 



ds > n> AT. 



Let us consider, for n > 1, r„ the stopping time 



r„ = inf | u > i : y Q< 

We get from the previous equation 

P t < E (^exp (J " e A ° (0c(0)6f(s) + ^ T (a{s) + 2/3 \Y^\)) ds^j P Tn | T t 

and, in view of the integrability assumption on a, Y and Y', since A s | < /3T, we can send n 
to infinity to obtain 



Pt < E f exp f y e As (0c(0)5/(s) + 7 e 0T (a(s) + 2/3 |Y S 'Q) ds j P T 
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On the other hand, (£ - 0f) = (1 - 0)f + (f - f) < (1 - 0)|f | + 6><5^, and since <5£ and <5/ are 
nonpositivc, wc finally derive the inequality 



exp 



1 <- 



f3T+A t 



(Yt - BYD 
y9 



5£ + / Sf(s) ds + 7 e 2 ^ T |f | + / (a(s) + 2/3 |Y S '|) ds 



1 -( 
< E I exp ■ 



In particular, since i5f and Sf are nonpositivc and f3T + A t > 0, 



(12) 



y - 017 < 



i - 



logE exp { 7 e 2/3T |f | + / (a(s) + 2/3 |Y S '|) ds 



and sending to 1, wc get Y t —Y^<0 which gives the first part of the result. 

For the second part of the theorem, if in addition we have Yq = Y ', then, coming back to 
the inequality (12), we deduce that 



exp ■ 



70 



1 - 



lei 



(a(s) + 2/3|Y s '|)ds 



8£ + J Sf(s) dsj + 7 e 2/3T 
and, once again, sending to 1, we get 

WW <fa=o ex P + ^ ("W + ^ l y s'l) d« 

which says that P (s£ = 0, J Q T <5/(s) ds = o) > 0. 



< e 7e " Ty ° < E 



□ 



As a byproduct, we obtain, using Corollary 4, the following existence and uniqueness result. 

Corollary 6. Let the assumption (A. 2) hold and let us assume moreover that the random 
variables |f| and \a\i have exponential moments of all order. Then the BSDE (1) has a unique 
solution (Y, Z) such that Y belongs to £ and Z belongs to M p for each p > 1. 



4. Stability. 



As we have seen in the previous section, the convexity of the generator with respect to the 
variable z leads to the comparison theorem. It turns out that we can also derive a stability 
result under this assumption. To be more precise, let us consider a generator / for which (A. 2) 
holds with parameters (a, /3, 7) together with a sequence of generators (f n )n>i '■ for each n > 1, 
f n satisfy (A. 2) with parameters (a„, /3, 7). Finally, let f and (£, n ) n >i be random terminal values 
such that, for each A > 0, 



E e MI?l+Mi) 



sup n>1 E 



M\U + \<* n \i) 



< 



(13) 



According to Corollary 6, let (Y, Z) be the solution to the BSDE (1) and let us introduce, 
for each n > 1, (Y™, Z n ) the solution to the BSDE 



yn 



f n ( S ,Y s n ,Z2)ds 



Z™-dB s . 
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Proposition 7. Let (A. 2) hold for f and f n and let us assume moreover that the inequality (13) 
holds true. 

If £n — ► £ P-o.s. anrf, m being the Lebesgue measure on [0,T], to ® P-a.e.,, /or eac/i 
(y,z) g R x R d , f n (t,y,z) — ► f(t,y,z), then, for each p > 1, 



E 



cxp (su Po < t < T \Y t n - Y t\) P +(f Q 



p/2 



ZXds 



Proof. It follows from Corollary 4 and the integrability assumptions (13) that the sequence 
((Y n , Z n )) n>1 satisfies, for eachp > 1, 



sup J 

n>l 



exp(su Po < t < T \Y t n \f + ^ 



P/2' 



\Z?\ 2 ds 



< +oo. 



It is thus enough to prove that 



su Po<t<r 



\Y t n ~Y t \+ f 
Jo 



\Zl l -ZA'ds 



converges to in probability to get the result of the proposition. 

Let us fix 9 £ (0, 1) and n > 1. First of all, since f n is convex in z and /3-Lipschitz in y, we 
can argue exactly as in the proof of the comparison theorem (Theorem 5) to estimate Y t n — 6Y t . 
Setting P t = e ceAtUt and Q t = cP t e A >{Z? - 6Z t ) where 



a t = [/„ (t, Yl\ Z t ") - /„ (t, 6Y t , Z?)] I [Y t n - 9Y t ] , A t = J a(s) ds, c = j—^ 

we have dPt = —Gt dt + Qt ■ dB t , with, see the inequality (11), 

G t < P t e At {c65 n f(t) + je^ T (a n (t) + 2(3 \Y t \)) , 

where 6 n f(t) stands for (/„ — /) (t, Y t , Z t ). Buy taking into account the fact that a(t) is bounded 
by (3, we have 

G t < Ptie^ T 
Let us introduce as before the process 

D? = exp ( 7 e 2/3T 

we get from Ito's formula 



\Snf(t)\ 



1 - 



+ a n (t) + 20\Y t \ 



{a n (s)+2[3\Y s \)ds 



D?P t < E D^P T 



ye 



2/3T i-T 



1 



P s D n s \5 n f{s)\ds 



Tt 



We observe that, since a is bounded by (3, 



Ps < G n (6) := sup exp 

0<t<T 



2/3T 



1 - 



\Yt\ + \Y t n \ 



Pt < X n {0) :=cxp 



7 e 



2/3T 



-dr-^|vie-6»e™i) 
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With these notations, we derive the inequality 



Y n -OYt < 



(1 - 9)e-^ T - A * 



7 e 



2/3T 



7 



logE \D^X n (6) + ^ ^D^G n {9) / \8 n f(s)\ds 



and finally, since log a; < a;, 
1 - 



Yt-Yt < (l-9)\Y t \ + ^—^E(D^X n (9)\T t ) + e 2l3T E \D%G n (6) jf \S n f(s)\ds | . (14) 

Now we want to find an upper bound for Y t — Y t n and to do this we use the same approach. 
Once again, let 9 6 (0, 1) and let us set U t = Y t - 9Y t n , V t = Z t - 9Z?. We have 

dU t = -H t dt + Vf dB t 

with H t = f(t, Y t , Z t ) — 9f n (t, Y™, Z"). To get the same inequality, we split H in the following 
way: 

H t = f(t, Y t , Z t ) - f n (t, Y t , Z t ) + f n (t, Y t , Z t ) - 9f n (t, Y t n , Z?) 

= -*„/(*) + f n (t, Y t , Z t ) - 9f n (t, Y t , Z") + 9f n (t, Y t ,Z?) ~ 9f n (t, Y t n , Z?). 



As before, the convexity of /„ with respect to z leads to 

f n (t, Y t , Z t ) - 9f„(t, Y t ,Z?) < (1 - 9)(a n (t) + 0\Y t \) + 

For the third term, let us introduce the process 



2(1 



/.\ fn(t, Y t , Z") - f n (t, Y t n , Z™) 

a(t) Y t -Y n i-m-Y^x) 

which is bounded by (3 so that 

9f n (t, Y t , Z?) - 9f n (t, Y t n , Z?) = 9a(t)(Y t - Y") = a(t.)(9Y t -Y t +Y t - 9Y t n ) 

< (3(1 - 9)\Y t \ + a(t)U t . 



It follows from the previous inequalities that 

H s < \S n f(s)\ + (1 - 9)(a n (s) + 2[3\Y S 



2(1 



^— -\V S \ 2 + a{s)U s 



with a bounded by (3. It follows from Ito's formula that, taking as usual A t = J ' a(s) ds, 



with 



e At U t = e AT U T + I e A °F s ds- I e A "V s dB s , 0<t<T, 



1 Wtf. 



F s < \6 n f(s)\ + (1 - 9){a n (s) + 2f3\Y s \) + 



2(1- 



This is the same inequality as (10). As a by product, we deduce that the inequality (14) is 
valid also for Y t — Y t n , namely, 



Y t -Y t n <(l-9)\Y t n \ 



1 



-E (D^X n (9) | T t ) + e 2/3T E \D^G n {9) J \5 n f(s)\ds 



Tt 
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and finally we have 

117" " Yt\ <(l-0) (\Y t n \ + \Y t \) + — E(fl?I„(e) | T t ) 



D^G n (9) / \S n f(s)\ds T t . (15) 



Let us fix e > 0. Wc have from (15) and Doob's maximal inequality 



sup \Y t n -Y t \ > e ) < ^ — ^ 

0<t<T J 



sup (|y t n | + |y t |) 

0<t<T 



3(1- 



7£ 



— E [D£X„(0)]+e 2 ^ 



■E 



£>fGn(0) / |<Jn/(s)|dfl 



and since, for 9 6 (0,1), the sequences (sup < t < T (|Y t n | + |Yt|)„>i), (£>t)„>i and ( G «(0)) n >i 
are bounded in all L p spaces, we deduce from Holder's inequality 



sup \Y t n -Y t \>e) < l -J-C{l + \\X n (9)h) + — ) 

0<t<T I 



\Snf(s)\ds 



(16) 



Let us recall that X n (9) = cxp ~ 0£| V |^ — #C l |)) so, as n goes to oo, X„(0) con- 

verges to cxp (7e 2 ^ T |£|) almost surely and actually in all L p spaces in view of the integrability 
assumptions on the sequence (£ n ) n >i (sec (13)). Moreover, \5 n f(s)\ = \f — f n \(s,Y s , Z s ) con- 
verges to to ® P-a.e. and, since 

\S n f{s)\ < a(s) + a n (s) + 20\Y,\ + 7 |^| 2 , 

Corollary 6 and the inequality (13) ensures the convergence to of \S n f(s)\ds in all L p 
spaces. The inequality (16) yields, for all < 9 < 1, 



n->+oo \0<t<T 



limsupP sup \Y t n -Y t \>e < C(l+||exp(7e 



2f3T\ 



and, sending 9 to 1, we obtain the convergence of sup 0<4<T |Y t " — Y t \ to in probability. As 
we mentioned before, we obtain the convergence of e sup °£ 4 £ T i( ~ Y * in all L p spaces. 

To get the convergence of Jq \Z™ — Z s \ 2 ds to in probability, let us just mention that, from 
Ito's formula wc have 



E 



\Z" - ZJ 2 ds 



< 



|e„-e| 2 + 2 sup \Y t n -Y t \ f ' \f n (s,Y s n ,Z")-f(s,Y s ,Z s )\ds 

0<t<T Jo 



from which the result follows directly. 



□ 



5. Application to quadratic PDEs. 



In this section, we give an application of our results concerning BSDEs to PDEs which are 
quadratic with respect to the gradient of the solution. More precisely, we want to obtain the 
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nonlinear Feynman-Kac formula in this framework. Let us consider the following scmilinear 
PDE 

dtu(t, x) + Cu(t, x) + / (t, x, u(t, x), a*V x u(t, xj) = 0, u(T,-)=g, (17) 
where £ is the infinitesimal generator of the diffusion X to - x ° solution to the SDE 



X t = x + f b{s,X s )ds+ f a{s,X s )dB s , t <t<T, X 

Jto Jt 



t = X , t < to. 



(18) 



The nonlinear Feynman-Kac consists in proving that the function defined by the formula 



V(t, x) G [0, T] x R ,; 



u(t, x) := y/ 



(19) 



where, for each (t ,x ) G [0,T] x R", (Y to < x ° , Z to ' x °) stands for the solution to the following 
BSDE 



Y t = g (X^ x °) + f(s, X***,Y„ Z s ) ds 



Q<t<T, 



(20) 



is a solution, at least a viscosity solution, to the PDE (17). Our objective is to derive this 
probabilistic representation for the solution to the PDE when the nonlinearity / is quadratic 
with respect to V x it and when g is an unbounded function. 

Let us first give our assumptions concerning the linear part of the PDE namely the coeffi- 
cients of the diffusion. 

Assumption (A.3). Let b : [0,T] x R" — > R™ and a : [0, T] x R" — ► R nxd be continuous 
functions and let us assume that there exists (3 > such that: 

(i) for all t E [0,T], |6(t,0)| < (3, and 

V(x,x') G R™ x R™, \b(t,x) -b{t,x')\ + \cr{t 7 x) -a(t,x')\ </3\x-x'\ ; 

(ii) a is bounded. 

Under the assumption (A.3), for each [to,Xo) G [0, T] x R", the SDE (18) has a unique 
solution denoted X to - x °. Classical results on SDEs show that, for each p > 1, X to ' x ° belongs to 
S p . Actually, since a is assumed to be a bounded function, for 1 < p < 2, we have 



VA > 0, 



E 



su Po<t<T e 



AX 



<CeW, 



where the constant C depends upon p, T, (3, A and ||<t||oo- Indeed, we have 



sup |X*°' X0 | < bo| +/3T + /3 / sup \X l t ' X0 \ds+ sup 



tn<t<u 



rtn.Xn 



to t <t<s 



tn<t<T 



a(s,Xl°> x °)dB s 



and we deduce from Gronwall's lemma the inequality 



sup \X! ' X0 \ < [ b 1 +(3T+ sup 



tn<t<u 



t <t<T 



f a(s,X t s °> Xa )dB s 

Jtn 



3 /3T 



It follows from the Dambis-Dubins-Schwarz representation of the continuous martingale 

/ a(s)dB s , with a(s) =a(s,X^ x °) 1* 0<s <t + W 
Jo 
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that, for A > and 1 < p < 2, 

sup e A|Bt|P , 
<lkllSo T 

which is a finite constant depending on p, A, T and ||<t||oo. 

With this observation in hands, we can give our assumptions on the nonlinear term of the 
PDE, the generator /, and the terminal condition. 

Assumption (A. 4). Let / : [0, T] x R" x R x R d — ► R and g : R" — ► R be continuous and 
let us assume moreover that there exist two constants f3 > and 1 < p < 2 such that: 

(i) for each (t,x,z) £ [0, T] x R" x R d , 

V(m')6R 2 , \f(t,x,y,z)-f(t,x,y',z)\<*3\y-y'\ ; 

(ii) for each (t, x, y) £ [0, T] x R" x R, z i — ► f(t, x, y, z) is convex on R d ; 

(iii) for each (t, x, y, z) £ [0, T] x R" x R x R d , 

\f(t, x, y,z)\ + \g(x)\ < d (1 + III" + |y| + |z| 2 ) . 

Since, for 1 < p < 2, sup l^* ' 310 | P has exponential moments of all orders, the growth 
condition on / and g allows one to use Corollary 6 to construct a unique solution, (Y tf > - x ° , Z tf > ' x ° ) , 
to the BSDE (20). Moreover, Y to ' xa belongs to £ and the process Z to ' x ° belongs to the space 
M p for cachp > 1. 

It is worth noticing that Y t °' x ° is actually deterministic for each point (to, %o) £ [0, T] x R". 
Indeed, the process {xl°' x °} t <t<T is known to be measurable with respect to the filtration 
generated by the increments of the Brownian motion after time to and it is by now well known 
that the process { (Y* 0,x ° , zl°' x °) } <t<T inherits this property. As a byproduct of this ob- 
servation, u defined by the formula (19) is a deterministic function. Let us prove that u is a 
viscosity solution to the PDE (17). 

Proposition 8. Let the assumption (A. 3) and (A. 4) hold. 

The function u defined by (19) is continuous on [0,T] x R" and satisfies 

V(t,a?) £ [0,T] xR", \u(t,x)\ <C(1 + |.t| p ). 
Moreover u is a viscosity solution to (17). 

Before proving this result, let us recall what is a viscosity solution to (17). 

Definition. A continuous function u on [0, T] x R™ such that u(T,x) = g(x) is said to be a 
viscosity subsolution (respectively supersolution) to (17) if 

d t ip(t ,x ) + £(p(t ,x ) + f {t , xo,u(t Q ,xo) 1 cr*V x ip(t ,x )) > 0, (respectively < 0) 

as soon asu-ip has a local maximum (respectively minimum) at (to, xq) £ (0, T) x R n where 
if is a smooth function. 

A viscosity solution is both a viscosity subsolution and a viscosity supersolution. 



E 



sup exp A 

t <t<T \ 



a (s, X 



to,Xo\ 



dB, 



< E 
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Proof. Our proof uses standard arguments since, in our setting, the BSDEs have all the prop- 
erties required: comparison and stability. First of all, the continuity of u is an immediate 
consequence of the stability property (Proposition 7) since the map (t, x) i — > X f > x is known to 
be continuous. Secondly, the growth of the function u comes directly from the general estimate 
on Y sec (5). 

Let us briefly explain why u is a viscosity subsolution to (17). Let p be a smooth function 
such that u — p has a local maximum at the point (to,xo) £ (0,T) x R™. Without loss of 
generality, we can assume that tp(to,xo) = u(to,Xo) and we want to prove that 

d t ip(t ,x ) + C(p(to,x ) + f (t ,xo,u(to,x ),a*S/ x (p(to,x )) > 0. 

Let us assume that the previous inequality does not hold. By continuity, there exist S > and 
< a < T - t such that 

w(f, x) < <p(t, x), et, dt<p(t, x) + C<p(t, x) + f (t, x, u(t, x), a* V x (p(t, x)) < — S 

as soon &s to < t < to + a and \x — a?o | < <x 
Let r be the following stopping time 

t = inf {u > t : {X^ '* - x | > a} A (t + a). 

The proof consists in applying the comparison theorem, actually the strict version, to the 
processes 

(Y%?\l t < T Z*>*°) and (<p (t A r, X%?) , l*< T cr*V^ (t, X*"*")) 
respectively solution to the BSDEs 

/to+a rta+a 
l s < T f{s,Xl^,Y s ,Z s )ds-J Z s 



dB 8 

Yl = <p (r, X***) + J° " -1 S < T {d t p + Cp} (s, X**°) ds - J° ° Z' s 



dB s 



In order to compare the terminal conditions and the generators of these BSDEs, let us recall 
that the uniqueness of solutions to (20) yields the Markov property: 

Y?°> x ° = u(t,X t t °' x °) , t>t . 

Hence, we can rewrite the first BSDE in the following way 

Y t =u(T,X t T °' x °)+ J° + 1 S < T f {s,Xi a ' x ° ,u{s,X^ Xa ) ,Z S ) ds- J° + Z s -dB s . 

By definition of the stopping time r, we have u (r, X* ^ ) < <p (r, X*°' x °) together with 
1 S < T f (s, Xl°> x °,u (s, Xl°' x °) , Z' s ) = l s < r / (s, Xl°' x °,u (s, Xt°> x °) , tr* V x¥ >(«, X^)) 

< -1,< T {d t <p + Lp] (s, X*»*<>) - 61 S < T . 

Moreover, it is worth noticing that 

-l.<r <Pt<P + Cp + f) (s, X^ x \u (s, Xl°> x °) , a*V x cp (s, X^ x °)) ds > 5{r - t Q ) 

and of course P(t = to) = 0. Thus, we can apply the second part of the comparison theorem, 
Theorem 5 above, and we get u{to,Xo) = Y to < Yl = tp(to,xo). But we have assumed that 
u(to,xo) = p>(to,xo)' u has to be a subsolution to (17). The fact that u is a supersolution and 
thus a solution can be shown in the same way. □ 
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Remark. When / does not depend on the variable y it can be shown that u is the unique viscosity 
solution with quadratic growth: \u(t, x)\ < C(l + |a;| 2 ). This follows from the uniqueness 
results in [5] concerning Bellman-Isaacs equation. 
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